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We report evidence for the decay B The analysis is based on a data sample of 152 

million BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB e"'"e~ storage 
ring. We detect a signal for — > n^n" with a significance of 3.4 standard deviations, and measure 
the branching fraction to be (1.7 ± 0.6(stat) ± 0.2(syst)) x 10"". 
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Recent measurements at B factories have improved 
significantly our knowledge of CP violation and heavy 
flavour physics. In particular, measurements of the 
mixing-induced CP violation parameter sin20i 0, Q 
lend strong support to the Kobayashi-Maskawa mecha- 
nism "s). It is of great importance to test precisely this 
theory with complementary measurements, such as those 
of the other unitarity triangle angles 02 and (ps 0]. 

The most promising technique for measuring (/>2 is by 
studyingtime dependent asymmetries in the i?° tt+tt^ 
system ]^ Q . The extraction of (l>2 from the observables 
is not trivial, however, because there are contributions 
from more than one weak phase. In order to disentan- 
gle 02, an isospin analyses of the tttt system can be per- 
formed 7]. One essential ingredient for these procedures, 
and hence for the measurement of 02, is knowledge of 
the branching fraction for the decay B — > ttOttO. Such 
knowledge would also play a pivotal role in the under- 
standing of charmless hadronic B decays. Previously, 
upper limits for the branching fraction of B ttOttO of 
(3 - 6) X 10"*^ have been reported [lIllIOl- Theoretical 
predictions are typically around or below 1 x 10~^ 
but phenomenological models incorporating large rescat- 
tering effects can accommodate larger values Il2j . The 
BaBar group recently measured this branching fraction 
to be (2.1 ± 0.6 ± 0.3) x lO'^ [H. 



In this paper, we report evidence for the decay B'^ — )■ 
TT°Tr°. The results are based on a 140 fb^^ (152 M BB 
pairs) collected with the Belle detector at the KEKB 
e+e" storage ring 0|. KEKB operates at a center-of- 
mass (CM) energy of y/s ~ 10.58 GeV, corresponding to 
the mass of the T(4S) resonance. Throughout this pa- 
per, neutral and charged B mesons are assumed to be 
produced in equal amounts at the T(4S), and the inclu- 
sion of charge conjugate modes is implied. 

The Belle detector is a large-solid-angle spectrometer 
consisting of a three-layer silicon vertex detector, a 50- 
layer central drift chamber (CDC), an array of threshold 
Cherenkov counters with silica aerogel radiators (ACC), 
time-of-flight scintillation counters, and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located 
outside of the coil is instrumented to detect mesons 
and to identify muons. A detailed description of the Belle 
detector can be found elsewhere 15]. 

Pairs of photons with invariant masses in the range 
115 MeV/c^ < nij^ < 152 MeV/c^ are used to form 7r° 
mesons; this corresponds to a window of ±2.5(t about 
the nominal tt*^ mass, where a denotes the experimental 
resolution, approximately 8 MeV/c^. The measured en- 
ergy of each photon in the laboratory frame is required 
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to be greater than 50 MeV in the barrel region, de- 
fined as 32° < < 128°, and greater than 100 MeV 
in the end-cap regions, defined as 17° < 8^ < 32° and 
128° < 0^ < 150°, where 6^ denotes the polar angle of 
the photon with respect to the beam line. 

Signal B candidates are formed from pairs of tt*' mesons 
and are identified by their beam energy constrained 
mass Mbc = -v/^blam^^P^ ^'^^ energy difference AE — 
Eg — i?boami where ii'boam denotes the beam energy and 
and Eg are the momentum and energy, respectively, 
of the reconstructed B meson, all evaluated in the e+e" 
CM frame. We require 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ 
and -0.3 GeV < A£; < 0.5 GeV. The signal effi- 
ciency of the kinematic reconstruction is estimated using 
GEANT-based ^ Monte Carlo (MC) simulations and 
found to be 26%. The resolution for signal is approxi- 
mately 4 MeV/c2 for Mbc and tloo MeV for AE. 

We consider background from other B decays and from 
e+e" ^ qq {q = u, d, s, c) continuum processes. A large 
generic MC sample shows that backgrounds from b ^ c 
decays are negligible. Among charmless B decays, the 
only significant background is 5+ p+Tr". We take 
these events, which populate the negative AE region, 
into account in the signal extraction described below. 

The dominant background is due to continuum pro- 
cesses. We discriminate signal events from the qq back- 
ground using the event topology. In order to increase the 
expected sensitivity to the signal, we have improved the 
continuum rejection technique used in our previous pub- 
lication 0. Previously, we have defined modified Fox- 
Wolfram moments [l2| that treat particles involved in 
the signal B candidate (s) separately from those in the 
rest of the event (o). We extend this idea, taking into 
account the missing momentum in the event, which we 
treat as a third category (m). We achieve some addi- 
tional discrimination by considering charged and neutral 
particles in the o category independently, and by taking 
the correlations of charges into account. In our previ- 
ously used continuum rejection technique, the moments 
are normalized relative to the zeroth moment; in this im- 
proved technique we do not normalize in this way. We 
combine 16 modified moments with the scalar sum of the 
transverse momentum into a Fisher discriminant |l8l |. 
and tune coefficients separately for seven categories of 
missing mass squared to maximize the separation be- 
tween signal and background. MC studies indicate that 
this redefinition of the Fisher discriminant leads to a 24% 
improvement in the maximum value of the figure of merit 
(FOM) defined as N^^'p/^/N^, where iVJ"^? and N°g§ 
denote the expected signal and observed background 
yields in a region 5.27 GeV/c^ < Afbc < 5.29 GeV/c^ 
and -0.20 GeV < AE < 0.05 GeV. 

The angle of the i?-meson fiight direction with respect 
to the beam axis (Ob) provides additional discrimination. 
A likelihood ratio TZg = Cs/{Cs + Cqq) is used as the dis- 
crimination variable, where Cs denotes the product of 



the individual Fisher and 9b likelihoods for the signal 
and £qq is that for the qq background. The likelihood 
functions are derived from MC for the signal and from 
events in the Mbc sideband region (5.20 GeV/c^ < Mbc < 
5.26 GeV/c^) for the qq background. We find additional 
discrimination between signal and background using the 
Belle standard algorithm for ^-flavour tagging [l', ^ . The 
fiavour tagging procedure yields two outputs: q — ±1 
(which we ignore), indicating the fiavour of the other 
purported B in the event, and r, which takes values be- 
tween and 1 and is a measure of the confidence that 
the q determination is correct. Events with a high value 
of r are considered well-tagged and are therefore unlikely 
to have originated from continuum processes. Moreover, 
we find that there is no strong correlation between r and 
any of the topological variables used above to separate 
signal from continuum. 

We combine r and TZs into a single multi-dimensional 
likchhood ratio (MDLR) defined as £f ^LRy^^MOLR _^ 
/:MDLR)^ ^l^gj.g £MDLR jenotes the likehhood deter- 
mined by the r-TZs two-dimensional distribution for the 
signal and C^^^ is that for the qq background. We 
then make a requirement on the likelihood ratio that 
maximizes the FOM. Incorporating the flavour tagging 
information in this way gives a 4% improvement in the 
FOM as compared to making a selection requirement on 
TZs alone. This criterion eliminates 99% of the qq back- 
ground while retaining 39% of the signal. The MDLR 
distributions for signal MC and for sideband data are 
shown in Fig. ^ We verified that there is no correlation 
between the MDLR cut and Affcc in the continuum MC 
events. 

The signal yield is extracted using an unbinned 
maximum-likelihood fit to the Mbc— Ai? two-dimensional 
distribution of the 596 candidates obtained after all the 
event selection requirements discussed above. The fit- 
ting function contains components for the signal, B~^ — s- 
p+7r° and qq background. The probability density func- 
tions (PDFs) for the signal and for B~^ p+Tr*' are 
taken from smoothed two-dimensional histograms ob- 
tained from large MC samples. For the signal PDF, dis- 
crepancies between the peak positions and resolutions in 
data and MC due to imperfect simulation of the tt" en- 
ergy are calibrated using — > tt^tt^ decay. The invari- 
ant mass distribution of Z?" are fitted with a bifurcated 
Gaussian for data and MC, and the observed discrep- 
ancies in peak position and width are converted to the 
differences of peak position and resolution of AE in our 
signal PDF, since observed differences are caused by im- 
perfect simulation of 7r° energy. For the D'^ daughter 
particles, similar momentum ranges and the same re- 
construction procedures are used as those for the sig- 
nal daughters. We find a 3 ± 9 MeV difference between 
MC and data for the AE peak position and a 35 ±12 
% discrepancy in the AE resolution. To obtain the two- 
dimensional PDF for the continuum background, we mul- 
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FIG. 1: Distributions of variables used for continuum sup- 
pression. (Top) Fisher discriminant using modified Fox- 
Wolfram moments for signal MC (solid line) and for events 
in the continuum dominated sideband region (dashed line); 
(bottom) MDLR for signal MC (solid line) and sideband data 
(dashed line). The dot-dashed line at MDLR = 0.95 indicates 
the selection requirement made. Apparent structure in the 
MDLR distribution arises from the several B-tagging modes 
that contribute to the r distribution. 
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FIG. 2: Result of the fit described in the text. (Top) 
Mhc projection for events which satisfy —0.2 GeV < AE < 
0.05 GeV; (bottom) AS projection for events which satisfy 
5.27 GeV/c^ < Mbc < 5.29 GeV/c^ The solid lines indi- 
cate the sum of all components, and the dashed, dotted and 
dot-dashed lines represent the contributions from signal, con- 
tinuum, and p'^vr". 



tiply the Ai? PDF, which is modeled with a first-order 
polynomial, with the Mbc PDF, for which we use the 
ARGUS function 19]. In the fit, the shapes of the signal 
and p+7r*^ PDFs are fixed, the normalization of 

£?+ p+tt'' is fixed according to the recent result from 
the BaBar collaboration , and all other fit parameters 
are allowed to float. The fit results are shown in Fig.|21 



The obtained signal yield is 25.6lg'4 with a statisti- 
cal significance (S) of 3.5, where S is defined as iS = 
^—21ti{Co/Cn^), and Cq and C denote the maximum 
likelihoods of the fits without and with the signal compo- 
nent, respectively. In order to take the uncertainty in the 
contribution from — > p+Tr*^ into account in the signif- 
icance calculation, we repeat the fit after increasing this 
contribution according to the error in the measurement 
of its branching fraction. In this case we find S — 3.4; we 
interpret this value as the significance of our signal. We 
also obtain consistent results when the normalization of 
the _B+ p+tt'^ component is allowed to float in the fit. 

We vary each calibration constant for the signal PDF 
by ±lcr and obtain systematic errors from the change in 



the signal yield. We also vary the 



P 



TT normal- 



ization, as described earlier, and assign systematic errors 
accordingly. Adding these errors in quadrature, we find 
the yield is 25.6l8 4(stat)l};4(syst). 

In order to obtain the branching fraction, we divide 
the signal yield by the reconstruction efficiency, measured 
from MC to be 9.9%, and by the number of BB pairs. 
The trigger efficiency is not corrected since it estimated 
to be 99% using our trigger simulator. We consider sys- 
tematic errors on the reconstruction efficiency due to pos- 
sible differences between data and MC. We assign a total 
error of 7% due to 7r° reconstruction efficiency, measured 
by comparing the ratio of the yields of the rj Tr^Tr^Tr" 
and r] —* 77 decays. The experimental errors on the 
branching fractions of these decays Q are included in 
this number. We compare the performance of the con- 
tinuum suppression requirement on a control sample of 
B+ D°{-^ ii:+7r"7r°)7r+ in data and MC; a systematic 
error of 2% is assigned. The efficiency of the MDLR cut 
for the MC control sample is close to that for the signal 
MC. Finally, we assign a systematic error of 0.5% due to 
the uncertainty in the number of BB pairs (152.0 ± 0.7) 
X 10^, and obtain a branching fraction of 



B{B° ^ 7r°7r°) = (1.7 ± 0.6(stat) ± 0.2(syst)) 



X 10 
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The result is stable for variations of the MDLR cut; for 
example, if we use 0.925 and 0.9 cut, we obtain 1.5 x 10^^ 
and 1.8 x 10~^, respectively, which are within the system- 
atic error. 

In conclusion, we have measured the branching frac- 
tion of B'^ 7r°7r° from a data sample of 152 mil- 
lion BB pairs collected at the T(4S) resonance with 
the Belle experiment. We obtain 25.6lg;4(stat)lJ;4(syst) 
signal events with a significance of 3.4 standard devi- 
ations (a). We measure the branching fraction to be 
(1.7± 0.6(stat) ± 0.2(syst)) x 10"^. This result super- 
sedes our previous result. We find evidence for B'^ 
tt'^tt'^ at a level above most theoretical predictions 0| 
and consistent with the BaBar measurement 0|. While 
larger values of the branching fraction for 7r°7r° 
enhance the feasibility of the isospin analyses 0, more 
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precise measurements of the branching fraction, in ad- 
dition to studies of the CP asymmetries in aU B — s- tttt 
decays, wih be required in order to determine 02 from 
the TTTT system. 
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